Impact of the strong electromagnetic field on the QCD eifective potential for 
homogeneous Abelian gluon field configurations 
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The one-loop quark contribution to the QCD effective potential for the homogeneous Abelian 
gluon field in the presence of external strong electromagnetic field is evaluated. The structure of 
extrema of the potential as a function of the angles between chromoelectric, chromomagnetic and 
electromagnetic fields is analyzed. In this setup, the electromagnetic field is considered as an external 
one while the gluon field represents domain structured nonperturbative gluon configurations related 
to the QCD vacuum in the confinement phase. Two particularly interesting gluon configurations, 
(anti-)self-dual and crossed orthogonal chromomagnetic and chromoelectric fields, are discussed 
specifically. Within this simplified framework it is shown that the strong electromagnetic fields 
can play a catalysing role for a deconfinement transition. At the qualitative level, the present 
consideration can be seen as a highly simplified study of an impact of the electromagnetic fields 
generated in relativistic heavy ion collisions on the strongly interacting hadronic matter. 
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I. INTRODUCTION 



The purpose of this paper is to study a poten- 
tial influence of the strong electromagnetic fields, 
eH ~ ^QCD' °^ ^^'^ QCD vacuum structure. Elec- 
tromagnetic fields with the strength of this order 
can emerge in relativistic heavy ion collisions. Be- 
fore proceeding, we have to decode our understand- 
ing of the stock phrase "QCD vacuum structure". 
In pure gluodynamics a physical vacuum can be 
characterized, first of all, by two invariants com- 
posed of gauge field: scalar gluon condensate {g^F"^) 
and pseudoscalar condensate {g^FF) (for instance 
see discussion in [l[). Since parity is not broken 
in strong interactions, the pseudoscalar condensate 
must be zero. Significance of the composite field 
g^FF becomes manifest in terms of topological sus- 
ceptibility. In QCD with quarks another conden- 
sate {mipil)) comes into consideration. Identifica- 
tion of gauge field configurations which are carri- 
ers of condensates and the method of their incorpo- 
ration into the formalism of quantum field theory 
can be seen as the most fundamental step towards 
understanding the mechanisms of confinement, chi- 
ral symmetry breaking and hadronization in QCD. 
This statement can be perceived as a kind of plati- 
tude since implicitly this step has to be assumed in 
all approaches dealing with configurations like center 
vortices, monopoles, instantons, etc. However, the 
feeling of having just a commonplace here relaxes 
if one identifies explicitly the point in the formalism 
where relevant condensates can be allowed or denied 
to be nonzero. As has been amphasizcd in a recent 



paper [2|, this point can be recognised in the choice 
of a functional space of the gauge fields to be inte- 
grated over in the QCD functional integral. In the 
Euclidean functional integral approach to quantiza- 
tion of the pure YM theory one starts with a symbol 



Z = N I DAcxp{^S[A]}, 



where the functional space J- of fields is subject to 
certain conditions, which can disable, in particular, 
the gluon condensate (requirement of finite classical 
action S[A\, for instance) or enable it and restrict 
the type of fields which can contribute to the con- 
densates. The character of fields in T has to be 
defined self-consistently on the basis of quantum ef- 
fective action. Enabling the gluon condensate means 
that gauge fields A't should satisfy 



T = {A 



lim — 

V->cx) V 



dWFUx)F",,{x)=B^}.{l) 



First of all, the requirement of nonzero condensate 
B^ ^ Q singles out fields S" with the strength which 
is constant almost everywhere in i?^, i.e. the part 
of i?^ where the field is inhomogenous has measure 
(4- volume) zero. The rest of deviations from homo- 
geneity can be treated as fluctuations in the back- 
ground of -BfJ. Separation of the long range modes 
i?" responsible for gluon condensate and the local 
fluctuations QJi in the background i?^, must be sup- 
plemented by the gauge fixing condition. The back- 
ground gauge condition D{B)Q = is the most nat- 
ural choice. At the formal level, the separation can 
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be achieved by the insertion of identity 
1 = / DB'S>[A, B] f DQ f DioSlA^ - Q'^ - B"^] 

x5[D{B'^)Q'^], 



Q n 



Al = Bl + Ql, 



(2) 



where Q are fluctuations of the gluon field with 
zero gluon condensate: Q G Q. Field B" are 
long range field configurations with, in general, the 
nonzero condensate: B € B. Performing the stan- 
dard Faddeev-Popov procedure one arrives at 

Z = N' DB DQdct[D{B)D{B + Q)] 

B Q 

xd[D{B)Q]exp{-S[B + Q]}. 

The character of long-range flelds has yet to be iden- 
tified by the dynamics of fluctuations Q. At the for- 
mal level, integral over Q defines an effective action 
for the long range part of the gluon field 

Z = N' f DBcxp{^S,b[B]}. 



Gluon fields B" , which correspond to the global min- 
ima of S'cft[B], dominate over the integral in the 
thermodynamic limit 1^ -> oo and define the phase 
structure of the system. First of all, one has to take 
a look at fields with just constant strength. There 
arc two different kinds of this type of fields: Abelian 
covariantly constant fields i?" = — ^n^B^^a;^ and 
non- Abelian constant vector potentials i?" = const. 
Unlike the former, non-Abelian fields are unstable 
against small perturbations Q° (for comprehensive 
discussion of the effective potential in pure Yang- 
Mills theory see[3,|j]. Pagels and Tomboulis studied 
an effective action for these fields within the con- 
text of scale anomaly @ , Woloshyn and Trottier 
attempted lattice calculation 0. All these calcula- 
tions indicated a minimum of the effective action at 
nonzero Abelian (anti-)self-dual field. Recently, the 
effective potential was calculated within the func- 
tional RG [3|- The result has also indicated a mini- 
mum of the effective action at the nonzero Abelian 
(anti-)self-dual field. In Q, the Landau-Ginsburg 
Lagrangian for pure Yang-Mills gauge fields invari- 
ant under the standard space-time and local gauge 
SU{3) transformations was considered. It has been 
demonstrated that for N^ = 3 a. set of twelve degen- 
erated minima of the action density exists as soon 
as a nonzero gluon condensate is postulated in the 
action. The minima are connected to each other by 
the Weyl group transformations associated with the 



color su{3) algebra and parity transformation. The 
presence of degenerated discrete minima in the La- 
grangian leads to the solutions of the effective equa- 
tions of motion in the form of the kink-like gauge 
field configurations interpolating between different 
minima. The homogeneous field with a kink defect is 
the simplest example of gluon configurations which 
arc homogeneous almost everywhere in R and sat- 
isfy the basic condition Eq. ^ . The spectrum of co- 
variant derivative squared D^ in the presence of the 
simplest solution, which interpolates between self- 
dual and anti-self-dual Abelian homogeneous fields, 
was estimated. This kink configuration can be seen 
as a domain wall defect separating the regions with 
self-dual and anti-self-dual Abelian gauge field. On 
the domain wall the gluon field is Abelian with or- 
thogonal to each other chromomagnetic and chro- 
moelectric fields. For the aims of the present study 
it is important that the spectrum of D^ orl^ in the 
(anti-)self-dual field is purely discrete with bound 
state type eigenfunctions while for the crossed or- 
thogonal fields the spectrum is continuous with the 
Landau level structure and the corresponding wave 
eigen functions. 

The eigenvalues and the square integrable eigen- 
functions of D^ for the (anti-)self-dual field are 



Kmkl{x) 

(j)a{x) 



4B (r + 1) 






1 



\/n\m\klllTr'^ 



where r = k + n for the self-dual field, r ~ I + n 
for the anti-sclf-dual field, /3^ and 7^ are related 
to a set of creation and annihilation operators (de- 
tails can be found in '^). The spectrum is discrete. 
In this background no color charged waves are en- 
abled, and there are no charged particle degrees of 
freedom. This is understood below as confinement 
of dynamical charged fields. 

Inside the infinitely thin domain wall placed at 
xi = with the chromomagnetic field directed along 
the y axis and the chromoelectric field along the z 
axis the charged scalar field displays a continuous 
spectrum similar to the Landau levels. The eigen 
functions square integrable over X3 take the form 

4'n{P2,Pi\x2,X:i,Xi) = Ciqi{~ipiXi ~ip2X2)Xn{Pi\xz), 

where the functions Xn are 

Xn{PA\x:i) = exp <^ -2V2B (^X3 + ^j 



The eigenvalues look like 

MpI.pI) = 2V2B{2n + 1) + pI + pI 

and correspond to the color charged quasiparticles 
with mass ni^ = 2^/2B{2n + 1) freely moving along 
the chromomagnetic field: 



pl=pl+ml. 



(3) 



The purely discrete spectrum and bound state 
(four-dimentional oscillator) eigen functions can be 
treated as confinement of color charged fields in the 
(anti-)self-dual homogeneous field (in the bulk of 
R^). Landau levels and wave eigenfunctions indi- 
cate the absence of confinement at the domain wall. 
In other words, charged particles arc localized at the 
wall. 

It should be noted here that accurate separation 
of the specific Abclian part B of general vector po- 
tential A as in Eq.(l2]) is a complicated problem. 
The methods to tackle the problem were studied in 

[iS[iMl. 

For completeness, we have to mention that in the 
context of center symmetry the dominance of the 
lumpy gauge field configurations was discussed in 
lattice calculations pi3l - l20{ . In paper [131, ^"^ effec- 
tive model of SU(2) gauge theory for the domain 
wall formation was considered. 

The idea of the dominance of the gluon fields 
which are (anti-)self-dual Abelian almost every- 
where turned out to be phenomenologically efficient. 
The model of confinement, chiral symmetry breaking 
and hadronization based on the ensemble of Abelian 
(anti-)self-dual fields was developed in a series of 
papers J2l| - J23f . In the model, the direction of the 
gauge field in space and color space, and the duality 
of the field are random parameters of the domains 
as well as pisitions of domain centers. All configu- 
rations of this type are summed up in the partition 
function. The domain model exhibits confinement of 
static (square law) and dynamical quarks (absence 
of poles in the propagators of color charged fields, 
discrete spectrum of the corresponding differential 
operator), spontaneous breaking of the flavour chi- 
ral symmetry, Ux{\) symmetry is broken due to the 
axial anomaly, strong CP violation is absent in the 
model. With a minimal set of parameters (meson 
masses, gauge coupling constant, gluon condensate 
and mean domain size) the model gives rather ac- 
curate results for meson masses from all different 
parts of the spectrum: light mesons including ex- 
cited states, heavy- light mesons, heavy quarkonia). 
The decay constants and some form factors were also 
calculated within the model. The above mentioned 
kink conflgurations have not been yet incorporated 
into the domain model directly but strongly moti- 
vate it. 




FIG. 1: Effective potential (in units of B /Stt ) as a 
function of angles 6 and ^ for the pure magnetic field 
H — .9B and (f> = x ■ The minimum is at S = 0, .^ = n /2. 



Strong electromagnetic flelds can emerge in rel- 
ativistic heavy ion collisions (24| - [26| . Interplay of 
strong electromagnetic fields and nonperturbative 
gluon fields are expected to be important for under- 
standing the dynamics of hadronic matter in heavy 
ion collisions. In particular, these fields can initi- 
ate such phenomena as chiral and vortic magnetic 
effects Hi, [13, [13. 

In this paper, we study an impact of the strong 
electromagnetic field on strong interactions in the 
context of lumpy or domain structured gluon fields. 
The one-loop quark contribution to the QCD effec- 
tive potential for the homogeneous Abelian gluon 
fields in the presence of homogeneous electromag- 
netic field is evaluated. Extrema of the potential as 
a function of angles between chromoelectric, chro- 
momagnetic and crossed orthogonal electromagnetic 
fields are analysed. In this setup the electromag- 
netic field is considered as an external one while 
the Abelian part of the gluon field represents do- 
main structured nonperturbative gluon configura- 
tions related to QCD in the confinement phase. It is 
shown that the quark contribution is minimal for the 
crossed chromoelectric and chromomagnetic fields 
orthogonal to each other, which can be treated as 
a catalyzing impact of strong electromagnetic fields 
on deconfinement in hadronic matter. It should be 
stressed that this result has a very indirect relation 
to the real physics of heavy ion collision since it does 
not take into account the temperature and density 
effect. The present extremely simplified calculation 
can play an instructive role for more realistic con- 



sideration. The main qualitative result of this paper could trigger a deconfinement transition in QCD. 
is an observation that strong electromagnetic field 

I 



II. ONE-LOOP QUARK CONTRIBUTION TO THE EFFECTIVE POTENTIAL IN THE 
PRESENCE OF ARBITRARY HOMOGENOUS ABELIAN FIELDS 



The one-loop contribution to the QCD effective potential is defined by the Gaussian integral over quark 
fields ^ with the covariant derivative which includes both the electromagnetic field and the Abelian homo- 
geneous gluon field 



g-yC/„ff(G) =_\f D^D^i,^! d^^^{^){^P-m)'4,{x) ^ 



C/cff(G) = -— In 



V det(?^— to) 

oo 

= — / d^xYx / dm! [5'(a;, ccJto') — iS'o(a;, a;|m')] , 



where S{x^y\m) is the fermion propagator in external gauge fields with mass m. The following notation is 
used 

Gij ~ £ijkHk, Gik = £k, 

n = nHgi + qU, £ = nEgi + 9E. 

Here electromagnetic fields arc denoted as H and E, and Hgi and Egi - chromomagnetic and chromoclectric 
fields of the same value Hgi = Egi — B, q \s quark electric charge. Trace includes sum of the elements of 
diagonal matrices n, m and q. Two invariants of the gauge fields are 

Tv. = —Gpi^i/G^i, = 'zyU + c ), Q~ G^^Gpi^y = tic. 

The quark propagator can be calculated analytically (see appendix |X| , 

S{x, y\m) = (m + ipx)H{x, y\m), 

with 

00 2 

H{x,y\m) = S{x - y) = e~^""^"-^---^ / ds . ^_ . ^^ - (4) 

m^+:p^ IGtt^ J smh(s^Qcr_) smh(s^Qcr+) 

P+ cosh(.|f - n\) + P_ cosh(s|£ + n\) \a,AG,. - G.^P^^^^^-^ \o,AG,. + g..] ''"^^^ t ^'^ 

^ \£ — '}i\ ^ \c ^'K\ 

{ y/Q^coth (.y^Q^) - y/Q^coth (.sy/Q^) 
^^^P| 4(a+-a_) ^""'^^ 

yo^coth (.yog - y^coth (.yg^ ^^^^ _ ^^^^^^ _ ^ 1 ^^^ 



1 1 Tl i I Q? 

-(i±75), CT^,. = — [7^,71^], CT± = — ii±yi-^ 



Using this propagator one gets for the effective potential 



C/,T(G)=Tr- 



Q f ds _„,2^ cosh(5|£ - H\) + cosh(5|g + ^1) 
sinh(s-^ Qcr^) sinh(s-^Qo'+) 



Here sq regularises the UV divergence of the integral, trace denotes summation over the elements of the 
diagonal color matrix n as well as quark charges q and masses m for all flavours under consideration. Using 
the identities 



cosh(s|£ ± "HI) = cosh(s(p+ ± p_)), \J Q.a± = p±, 



(6) 



^/2 
one arrives at the renormalized effective potential written in terms of invariants TZ and Q 



U^siG) = Ul^^{G) + 5C/eff, <5C/eff 



B^ f ds 



C^eT(G) = ^ I -Tr„ 



87r2 



Tr-n 



i J J s 

So 



ds 



sx-^ coth(s>i:+)s>f_ coth(s>f_) — 1 5"(^+ + ^- 

o 



(7) 



K± 



V2B 



Vn + Q±Vn~Q 



Thr oug h the identity ([6]) this expression can be reduced to the well known form of the effective potential, 
see [ll[ and references therein. 

For the experimantal situation of heavy ion collisions, which we bear in mind, it is sufficient to consider 
the electric and magnetic fields orthogonal to each other and choose coordinate system with the z-axis along 
the magnetic field and the x-axis along the electric field, 

H.^SaH, E.^SaE, EH - 0, 

For the physical electric field orthogonal to the magnetic field invariants TZ and Q read 

n = {H^ -E^)/2 + n^B^ +hBHcos{0)+iBEcos{x)sin{^), 



Q = hBH cos{^) + ifiBE sm{e) cos(0) + n^B^(sin(6') sm{£_) cos {(j) - x) + cos{0) cos {£_)), 



(8) 



r 



where (0, 9) are the spherical angles of the chromo- 
magnetic field, and (x, ^^^ the spherical angles 
of the chromoelectric field in the chosen coordinate 
system. 



III. MINIMA OF THE EFFECTIVE 
POTENTIAL 

For the pure magnetic field {E = 0) the renor- 
malised effective potential is real. It depends on 9, 
^ and the difference between (/) and x- Figure [1] rep- 
resents the effective potential for the case of three 
quark flavours, the masses are taken the same for all 
flavours. A minimum is achieved at = 0, ^ = 7r/2, 
i.e. for the crossed orthogonal chromomagnetic and 



chromoelectric fields. The value of the potential does 
not depend on (/) and x- ^^ other words, the chromo- 
magnetic field is collinear to the magnetic field, and 
the chromoelectric field is orthogonal to the chro- 
momagnetic field. The polar angle x of Egi is not 
fixed. 

If the electric field is nonzero, the effective poten- 
tial becomes comlex. The analytical properties of 
the effective potential as a function of complex Q 
and TZ can be studied by means of the convergent 
series representation of the integral ([7]) obtained in 
|ll|. The nonzero imaginary part of the potential 
would mean instability in the system. However, the 
imaginary part does vanish at the pont in the space 
of angles where the real part is minimal. This can be 
verified by inspection through the straightforward 



calculation of the real and imaginary parts of the 
potential. The simplest way to find this stable min- 
imum of the potential is to require that the RHS 
of Eqs.® for TZ and Q to be real, which can be 
achieved by restriction imposed on the angles 

sin(^) cos(x) = 0, 
sin(6') cos{(j)) = 0. 

Figures [5] and |3] illustrate that a globall minimum is 
achieved at 6* = 0, ^ = 71-/2, x = 7r/2 or 3tt/2. It 
does not depend on (j). This minimum corresponds 
to the crossed orthogonal to each other chromoclcc- 
tric and chromomagnetic fields, the chromomagnetic 
field is collinear to the magnetic field, and the chro- 
moelectric field is orthogonal to both magnetic and 
electric fields. The minimum exists for any value of 
the electromagnetic fields. A value of the potential 
at the minimum depends on their strengh. 



IV. DISCUSSION. 

We have analysed dependence of the one-loop 
quark contribution to the QCD effective action on 
the angles between physical homogeneous electro- 
magnetic field and Abelian homogeneous Euclidean 
gluon field. The main result consists in observa- 
tion that in the presence of the external orthogo- 
nal magnetic and electric fields there exists a global 
mimium corresponding to the crossed orthogonal 
chromomagnetic and chromoelectric fields. Unlike 
the (anti-)self-dual homogeneous gluon field min- 
imisingthe effective potential in pure gluodynam- 
ics [H, 13, S Hi th^ crossed orthogonal gluon field 
does not support confinement of quarks: the color 
charged quasi-particles do exist and can move along 
the direction of the magnetic field, see Eq.Q . In 
this sense, strong electromagnetic field can rearrange 
the structure of the global minima of the effective 
action of QCD. 

If to take a liberty to extrapolate this extremely 
simplified setup to the situation of relativistic heavy 
ion collision, then one can expect that a strong 
electromagnetic field generated during the collision 
triggers quark dcconfinemcnt transition in hadronic 
matter. Within the context of the domain model 
and the kink solution one can think that a strong 
electromagnetic field produces a domain wall defect 
in the confining gluon background exactly in the re- 
gion where collision occurs. If so then deconfined 
quarks will move preferably along the direction of 
magnetic field but this will happen due to the gluon 
field configuration, due to QCD interaction, even af- 
ter the switching the electromagnetic field off. How- 
ever, prior to discussing the phenomenologically rel- 
evant observables arising from this effect, we should 
make the computation more realistic. First of all. 



eff 




YJ ren 
eff 




FIG. 2: Effective potential (in units of B^/Stt^) for the 
electric E = .55 and the magnetic H — .9B fields as 
a function of the angles 9 and C (0 — X = '"'/S - upper 
figure, cj) = 7r/2, x = 37r/2 - lower figure ). 



this would imply incorporation of the temperature 
and baryonic density effects, and taking an inho- 
mogeneity of the gauge fields into account as far as 
possible. 
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one can get 
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FIG. 3: Effective potential (in units of B^/Stt^) for the 
electric E = .5B and the magnetic H = .9B fields as a 
function of the angle 9 and ^ — -k 12 {(t> = X ~ ^/2 )• 



Appendix A: Quark propagator in the presence 
of arbitrary homogenous Abelian fields 

The quark propagator 

S{x, y\m) = ^l^^i^ - 2/) = ("^ + ipx)H{x, y\m) 



1 ~ smh{s\£ + n\) 



The integrand in Eq. (j4]) can be computed by 
means of the path integral representation, 

e^p{sD'^}d{x -y)= 6a^Pp exp i - / dl3a^{l3) 



+2^s / dfia^{fi)DA5{x~y) 



1 



IGtt^s^ 



exp 



[x ~ yf 

As 



l^x^G^^y^, I Y[ Vdet 0(fc) 



fc=i 



2s 



*^^P 1 X! r^Gp-ix - y)Oj{n)Gr,{x - y) ^ , 



H{x,y\m) 



1 



;S{x-y), 



can be represented in the form 

oo 

H{x,y) = I ds e"™'" e'^'^"^^ e'°^ 5{x - y)iKl) 



with 



where 

0„a{n) = 
and the determinant is 
detO(fc) =det 



^^ ai/^ na 






2\c2 



1 



7r4fc4 + 



7r2fc2 



With the notation 



s = (fH), 7^ = -(£2 + H2), 



-± = | i±V'-|^ ' -+^^- = 1' 



Using relations 



one arrives at 



[gG) ^{£- nr'p+ + {£ + nr'P- , 



det 0{k) 



Q.sV_ 



7r2fc2 



7r2/fc2 



(ctG)P± = -cr^„[G^,y =F G^ix], 



[] det 0^/2 (^) 



Lfc=0 



s2Q 



sinh(s ■^Qcr_) sinh(s-^Q(T-|_) 



-P± — 2''"'" ^'^'s)' ^M" ^ ■::£ p,vafiG ap 



The sum in the exponent with O J can be evaluated as 



2s 



J2 -^^G^{x - y)0-^{n)G,{x - y) = G^{x - y)G^{x ~ y)Ei + {x - yfll^, 



where the terms Si and S2 are given by 



E 



2sTT^n 



2^2 



2cr_s 



E 



2(7+ s 



E 



—1 ^ ' n— 1 ^ ^ n— 1 



li 






sQ 



E 



2 ^^ TT^n^ + 2s27^2r^27^ + Q2s4 2(cr+ - cr_) ^-^ 

n— 1 ^ ^ n— 1 

y^Scr+coth (s^Qcr_) - ^Qcr_ coth {sy/Qa+) 1 



7r2„2 + siQ ^2„2 + £ie 



<^+ 



4(cr+-cr-) 

as a result of this calculation one gets representation 

1 Q 



4s 



exp{sZ? }S{x — y) = 



167r2 sinh(s ^Qo-_) sinh(s^Q(T+) 



cxpi -T^^aGa^y^ 



X cxp ■ 



^Qcr+ coth (s^Qcr_) - ^Qcr_ coth (s^Q^ 



4(o-+-o-_) 



-(2^- -2/)' 



VQcr+COth (s^Qcr+) - ^Qcr_COth (s^Qcr_) 

4Q(f7+-(7_) 



G^i.G^p(a;-y)i.(a:-y)f 



Substitution of Eqs. (jA2p and (jA3P to Eq. (|A1|) leads to the quark propagator ^. 

I 



(A3) 



[1] P. Minkowski, Nucl. Phys. B177 (1981) 203. 

[2] L. D. Faddeev, [arXiv:0911.1013 [matli-ph]]. 

[3] P. Minkowski, Phys. Lett. B 76 (1978) 439. 

[4] H. Leutwyler, Phys. Lett. B96 (1980) 154; ihid 

Nucl. Phys. B 179 (1981) 129. 
[5] H. Pagels, and E. Tomboulis, Nucl. Phys. B 143 

(1978) 485. 
[6] H. D. Trottier and R. M. Woloshyn, Phys. Rev. Lett. 

70 (1993) 2053. 
[7] A. Eichhorn, H. Gies and J. M. Pawlowski, Phys. 

Rev. D 83, 045014 (2011) [Erratum-ibid. D 83, 

069903 (2011)] [arXiv:1010.2153 [hep-ph]]. 
[8] B.V. Gahlo and S.N. Nedelko, Phys. Part. Nucl. 

Lett., 8 (2011) 67 [arXiv:hep-ph/1006.0248v2]. 
[9] Y.M. Cho Phys. Rev. Lett. 44 (1980) 1115. 
[10] Y. M. Cho, J. H. Kim and D. G. Pak, Mod. Phys. 

Lett. A 21 (2006) 2789. 
[11] Y.M. Cho, and D.G. Pak, Phys. Rev. Lett. 86 

(2001) 1947. 
[12] S.V. Shabanov, J. Math. Phys. 43 (2002) 4127 [hep- 

th/0202146]. 
[13] S. V. Shabanov, Phys. Rept. 326 (2000) 1 

[arXiv:hep-th/0002043]; S. V. Shabanov and 



J. R. Klauder, Phys. Lett. B 456 (1999) 38. L. V. 

Prokhorov Yad. Fiz., 35 (1982) 229. 
[14] L. D. Faddeev, A. J. Niemi, Nucl. Phys. B 776 

(2007); ihid, Phys. Lett. B 449 (1999) 214. 
[15] Kei-Ichi Kondo, Torn Shinohara, Takeharu Mu- 
rakami, Prog. Theor. Phys. 120 (2008) 1 

[arXiv:0803.0176 [hep-th]]. 
[16] P. J. Moran and D. B. Leinweber arXiv:0805.4246 

[hep-lat] . 
[17] P. J. Moran and D. B. Leinweber PoS LAT2GG7 

(2007) 383 [arXiv:0710.2380 [hep-lat]]. 
[18] Ph. de Forcrand, A. Kurkela and A. Vuorinen Phys. 

Rev. D 77 (2008) 125014. 
[19] P. de Forcrand AIP Conf. Proc. 892 (2007) 29 

[arXiv:hep-lat/0611034]. 
[20] E. M. Ilgenfritz, et al, PoS LAT2007 (2007) 311 

[arXiv:0710.2607 [hep-lat]]. 
[21] G.V. Efimov, and S.N. Nedelko, Phys. Rev. D 51 

(1995) 176; Ya. V. Burdanov, et. al, Phys. Rev. D 

54 (1996) 4483. 
[22] A.C. Kalloniatis and S.N. Nedelko, Phys. Rev. D 64 

(2001) 114025; 
[23] A.C. Kalloniatis and S.N. Nedelko, Phys. Rev. D 66 



(2002) 074020; ibid, Phys. Rev. D 69 (2004) 074029; 

Erratum-ibid. Phys. Rev. D 70 (2004) 119903; ibid, 

Phys. Rev. D 71 (2005) 054002; ibid, Phys. Rev. D 

73 (2006) 034006. 
[24] D. E. Kharzeev, L. D. McLerran, and H. J. War- 

ringa, Nucl. Phys. A 803 227 (2008) 
[25] V. Voronyuk, V. D. Toneev, W. Gassing, 

E. L. Bratkovskaya, V. P. Konchakovski and 

S. A. Voloshm, arXiv:1103.4239 [nucl-th]. 



[26] V. Skokov, A. Y. lUarionov and V. Toneev, Int. 
J. Mod. Phys. A 24 (2009) 5925 [arXiv:0907.1396 
[nucl-th]]. 

[27] D. E. Kharzeev and D. T. Son, Phys. Rev. Lett. 
106, 062301 (2011) [arXiv: 1010.0038 [hep-ph]]. 

[28] O. Rogachevsky, A. Sorin and O. Teryaev, Phys. 
Rev. C 82, 054910 (2010) [arXiv: 1006. 1331 [hep- 
ph]]. 



